Functional Renormalization for Chiral and J7a(1) Symmetries at Finite Temperature 
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We investigated the chiral symmetry and (7a(1) anomaly at finite temperature by applying the 
functional renormalization group to the SU(3) linear sigma model. Expanding the local potential 
around the classical fields, we derived the flow equations for the renormalization parameters. In 
chiral limit, the flow equation for the chiral condensate is decoupled from the others and can be 
analytically solved. The Goldstone theorem is guaranteed in vacuum and at finite temperature, 
and the two phase transitions for the chiral and Ua(1) symmetry restoration happen at the same 
critical temperature. In general case with explicit chiral symmetry breaking, the two symmetries 
are partially and slowly restored, and the scalar and pseudoscalar meson masses are controlled by 
the restoration in the limit of high temperature. 



1. INTRODUCTION 

The chiral symmetry is one of the fundamental proper- 
ties of quantum chromodynamics (QCD). In chiral limit 
with zero current quark mass, the QCD Lagrangian den- 
sity respects the symmetry of Z7l(3) x Ur(3) = C/y(l) x 
Ua(1) x SUy(3) x SUa(3) at classical level. In vac- 
uum, the condensate of right-hand quark and left-hand 
anti-quark ((JlQr) breaks the SUa(3) symmetry, and the 
[7^(1) is broken by the anomaly due to the nontrivial 
topology of principal bundle of gauge fieldQ. The chi- 
ral symmetry breaking leads to a rich meson and baryon 
spectrum, and supplement the 17a (1) anomaly ex- 
plains the non-degeneracy of 77 and 77' mesons @. As 
a strong interacting system should approach its classic 
limit at high temperature, the chiral symmetry is be- 
lieved to be restored in hot medium. The relation be- 
tween the two symmetries in vacuum and at finite tem- 
perature has been studied for a long time. While the 
two symmetries are expected to be restored Q at high 
temperature, it is still an open question whether the two 
phase transitions happen under the same condition. 

The lattice simulation is a powerful tool to study the 
QCD symmetries in vacuum and at finite temperature. 
The chiral condensate is observed to decrease with in- 
creasing temperature and the chiral susceptibility shows 
a peak at a critical temperature T c . By a proper defi- 
nition, the topological charge and its susceptibility are 
used to describe the 17^(1) anomaly in the pure gauge 
field theory and the unquenched theory [H, @. In both 
cases, the susceptibility drops down above the critical 
temperature T c of the chiral restoration and the charge 
keeps an obvious deviation from zero at high tempera- 
ture T > T c . The simulation for the instanton model 
shows such a partial restoration too@. On the other 
hand, the observation of the hadron spectra provides an 
experimental way to test the restoration of the two sym- 
metries in hot medium created in relativistic heavy ion 
collisions. The mass shift due to the chiral restoration 
enhances or reduces the hadron thermal production, for 
instance for the kaon yields and ratios pHS!- The partial 
restoration of the [7a (1) symmetry is closely related to 
the production of 77 meson and spin-excited hadrons in 
hot medium [13, [H| ■ 



The effective modelsj!2l 1 1 31] are often used to describe 
the QCD phase structure, especially at finite baryon den- 
sity where the lattice simulation meets the fermion sign 
problem and can not yet present precise result. For 
the study of chiral symmetry at low energy, the effec- 
tive models at hadron level can include only scalar and 
pseudoscalar mesons. In this work we choose the lin- 
ear sigma model which has been widely discussed in 



vacuum|14j|. The simplest version of this model is with 
the symmetry of SU L {2) x SU R (2) or 0(4), including 
only the Goldstone modes 7r and their partner, the a 
meson. The linear sigma model is considered as a good 
laboratory for various approximation methods like mean 
field and the Cornwall- Jackiw-Tomboulis (CJT) approx- 
imation. In the frame of this model, the a and 7r prop- 
erties have been deeply investigated in vacuum and hot 
medium [l^, [l^]. The model is also extended to include 
fundamental represented quarks coupled to the mesons in 
Yukawa form. However, in the SU(2) version the 17a (1) 
anomaly can not be properly studied, since there are no 
enough flavors for r\ and 77' mesons which are related to 
the anomaly. In order to study the chiral symmetry and 
17a(1) symmetry at the same time, the flavor symmetry 
group is chosen as SUl(3) x SUr(3) in this work. 

The pion, kaon and eta mesons are the Goldstone 
modes corresponding to the spontaneous SUa(3) sym- 
metry breaking. These zero mass mesons are guaranteed 
by the Nambu-Goldstone theorem in vacuum as well as 
at low temperature T < T c . However, in the mean filed 
approximation[l5j and the CJT approximation[l|| Il7j |. 
the Nambu-Goldstone theorem is destroyed at low tem- 
perature. In this work we use the functional renormaliza- 
tion group (FRG) method to study the chiral symmetry 
and 17a(1) symmetry in the SU(3) linear sigma model at 
finite temperature. As a non-perturbative method [l8l- 
l20j , the FRG has been used to study phase transitions in 
various systems like cold atom gas[2l[, nucleon gas[22j 
and hadron gas p3l - [28j . By solving the flow equation 
which connects physics at different momentum scales, the 
FRG shows a great power to describe the phase transi- 
tions and the corresponding critical phenomena, which 
are normally difficult to be controlled in the mean field 
approximation because of the absence of the quantum 
fluctuations. Instead of adding hot loops to the thermo- 
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dynamic potential in usual ways of going beyond mean 
field, the FRG effective potential at tree level of mean 
field approximation includes already quantum fluctua- 
tions through the mass and coupling constant renormal- 
ization and can guarantee the Nambu-Goldstone theorem 
in the chiral symmetry breaking phase. 

We proceed as follows. In Section [31 we briefly re- 
view the SU (3) linear sigma model, following the nota- 
tion in [ItJ , and apply the functional rcnormalization to 
the model. The flow equations are derived in general 
case with explicit chiral breaking terms, and the one for 
the chiral condensate can be solved analytically in chiral 
limit. In Section [3] we present the numerical results for 
the light and strange quark condensates and the topologi- 
cal susceptibility at finite temperature. Finally in Section 
Uwc summarize our results. 



2. APPLICATION OF FUNCTIONAL 
RENORMALIZATION TO THE SU(3) LINEAR 
SIGMA MODEL 

The SUl(3) x SUr(3) linear sigma model has been 
widely studied in mean field and CJT approximation in 
vacuum and at finite temperature [13, [H,!!^]. Following 
the notations in [Ij} , the Lagrangian density of the model 
is expressed as 

C = Tr (fy^fl'V - mV>) +c[Det(<£)+Det(<£ t )] 



Ai [Tr 



A 2 Tr(^0) 2 + Tr[ff(0 + 0t)] , 



(1) 



where the meson matrix = T a <p a and the trace Tr are 
defined in the flavor space, the meson field 4> a = a a + in a 
contains the scalar part a a and the pseudoscalar part 7r Q , 
the 3x3 Gell-Mann matrices T a = A a /2 for a = 1, ■ ■ ■ ,8 
and T = ^/2/3 for a = obey the relations Tr(T a T 6 ) = 
5 ab /2, [T a ,T b ] = ifabcT c and {T a ,T 6 } = d abc T c with the 
structure constants f abc and d abc , m 2 is the mass param- 
eter, and c, Ai and A2 are the coupling constants. 

The Lagrangian density (JTJ) is invariant under the 
SUl(3) x SUr(3) transformation, except the last term 
which explicitly breaks the chiral symmetry, 



Tr[ff(0 + </)+)] = h a a a , 



(2) 



where the matrix H is defined as H = h a T a with 9 pa- 
rameters h a . 

The determinant term in ([I]) explicitly breaks the 
Ua(1) symmetry which in QCD is violated by the 
anomaly. If the coefficient c of the J7a(1) anomaly term 
vanishes, the symmetry group of the system is enlarged 
to U L (3) x Ur(3). 

In vacuum and at finite temperature but zero density, 
there are only scalar condensates 



{<!>)= T a {a) 



(3) 



where (X) means the ensemble average of the operator 
X. To simplify the notation, we use a to replace (a) in 



the following. Making a shift for the meson field <f> — s- 
(4>) + 84> and substituting it into the Lagrangian density 
(JTJ, the effective potential of the system at classical level 
can be written as 



m -2 n - - - 1 
U{(7) = -^-CTq — Gabc&a&bPc + ^Fabcd&aCrbCrcCrd 

— h„a n , 



(4) 



and the dynamical masses generated by the condensates 
can be extracted from the coefficients of the term (S(p) 2 
and form two 9x9 matrices Mg and Mp for the scalar 
and pseudoscalar mesons, 



(MlU 



(M 2 P ) ab 

with the coefficients defined as 



m 2 S a b — 6G a bcCTc + 4F abcd a c a d , 
m 2 S ab + 6G abc a c + 4H abcd a c a d 



(5) 



c r 9 

Gabc — 7 d abc + -dooo8ao8 b o8 c o 

L 2 
3 

— — (Saod-Obc + Sbod a o c + 8 c od ab o) 

F abcd = -f(S ab S cd + S ad S bc + 5 ac S bd ) 
A 2 

~t~ "7T" (,d a bed ecd ~t- d ade d eb c ~\~ d ace d ebd *) , 
o 

Habcd = -^SabScd 
A 2 

+ fadedebc + facedebd) ■ (6) 

o 

The physical condensates are determined by minimiz- 
ing the potential, 



dU(a) 







which leads to the gap equations 
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m 2 a n 



3G abc a b a c + -F abcd a b a c a d — h a 



0. 



(7) 



(8) 



In different thermodynamic environments, the conden- 
sate (4>) can be further simplified. Since the isospin zero 
mesons Co and 03 or and ay can couple to the vac- 
uum without violating Lorentz invariance and parity, the 
classical field matrix (</>) and the coefficient matrix H in 
the chiral breaking term contain only two components, 
((f)) = T a + T 8 a s and H = T h Q + T s h 8 . In order to 
simplify the expressions, we make a rotation for the con- 
densates oo and as and the chiral breaking parameters 
ho and hg, 



h a 



^0 

CTO 

h 
ho 



0-8, 

h 8 , 
h 8 . 



(9) 
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In terms of the rotated condensates a u and er s and the 
rotated breaking parameters h. u and h s , the classical po- 
tential and the gap equations can be explicitly expressed 
as 



TTI- - \ " l 1-2 

U{cr u ,a s ) = — {a u 



2V2 



^(al + ^) 2 + ^(at +2at), (10) 



and 



m 2 a v 



m 2 a s 



1 
1 

2V2 



ca u a s + Ai (a 2 , + a 2 ) a u + -X 2 crf. = K 



<*i + Ai (o* + a 2 ) a s + X 2 a 3 s = h s . (11) 



For the meson mass matrices Ms and Mp, each has only 

r2 - M« 2 n and 
r2 "JI/T2 



one independent off-diagonal element Mq S — iu 80 



four independent diagonal elements Mq , Mg 8 and 





= (Ml)u 


= (Af|) 22 


= (M|) 3 3, 


(12) 


m l 


= (M 2 ) 44 = 


= (Ml) 68 = 


= (^1)66 = 


(M|) r7 , 


m 2 


= (m£) u 


= (M P ) 22 


- (M P ) 33 , 




m 2 K 


= (M P ) 44 


= (M P ) 55 


- (Mp) 66 = 


= (M 2 P ) 77 , 



and diagonalizing the meson subspace a = 0, 8 generates 
the pseudoscalar mesons 77 and 77' and the corresponding 
scalar mesons. There are six parameters in the model, 
the mass m, the three coupling constants c, Ai and A 2 
and the two chiral breaking parameters h u and h s . They 
should be determined by the experimental data in vac- 
uum. Firstly, the partial conservation of axial- vector cur- 
rent (PCAC) leads to a relation between the condensates 
and the pion and kaon decay constants f T and fx, 



-U + Vk 

V2 ' 



(13) 



then the gap equations (jlip can be reexpressed by the 
Goldstone modes ir and K, 



h u = ml fir, 
, -m'lU + 5m 2 K f K 

V2 



(14) 



and the combination of the isospin zero pseudoscalar 
mesons 77 and 77' determines the couplings A 2 and c, 



to, c, Ai, A 2 , h u and h s are fitted by the experimental 
values of , m^ , f n , , and m 2 + m 2 , and one of the 
scalar meson masses. It is necessary to note that, in this 
way one can fix a group of parameters, but the obtained 
condensates may not correspond to the minimum of the 
potential. One should check the secondary derivative of 
the potential, d 2 U jdoidoj > 0. 

We now apply the functional renormalization group to 
the SU(3) linear sigma model. The core quantity in the 
frame of FRG is the averaged effective action at a 
momentum scale k in Euclidean space, 

r k [(<t>}}= [ d*x[Tr(Z k d li {<j>)1&'(<i>))+U k {(<i>)) +•■■], 



(16) 

where Z k is the wave function renormalization constant, 
Uk((4>)) is the classical potential (|TU1) but with renor- 
malized mass and coupling parameters mk,Ck,\\k and 
A 2 fc and scale dependent condensates a uk and a skl and 
the symbol ■ ■ • stands for the high order terms of the 
field ((f)) . The scale dependence of the averaged action is 
characterized by the flow equation [l8H2Q[ in momentum 
representation, 



dk 



p 



Tr 



(2) 



- 1 dR k 
dk 



(17) 

(2) 

where 1^ is the second order functional derivative of the 

averag cd action rj; 2) [(^}] = 5 2 T k / 5 (4>) 2 , and the infrared 
cutoff function R k which is used to suppress the quan- 
tum fluctuations at low momentum p < k is chosen as 
the optimized regulator function R k = (k 2 —p 2 )6(k 2 —p 2 ) 



3 lj | . From our numerical calculation shown in the next 
section, the symmetry restoration and meson mass spec- 
tra at finite temperature are not sensitive to the choice 
of the cutoff function (32|. 

Following the effective action flow starting from the ul- 
traviolet momentum k = A, the physics we are interested 
in could be obtained at k = 0. 

By assuming the space-time independence of the clas- 
sical field (</>), the effective action to the lowest order is 
determined by the classical potential only, 



d A xU k ((0)) 



(18) 



With the known infrared cutoff function R k , after doing 
the three-momentum integration for the mesons at finite 
temperature, the FRG flow equation can be simplified as 



A, 



| \/6a s m 2 K - \/2a ml + a s (m 2 
' (^+4ct s >8 

- - X 2 (2f K - U). 



Ik - A 



(15) 



d k U k 



67T 2 

6tt 2 



1 



Vsk \ 
6 ) 



r^TrD Sfe 



(19) 



Substituting the above obtained parameters into any of 
the two gap equations, one can get the relation between 
the mass to and the coupling Ai. By fitting a scalar 
meson mass, for instance m a , one can then separately 
fix the two parameters. In summary, the six parameters 



with the meson propagators D 



and D p ], 



J p k 



-kd k Z k 1 /Z k , where 



Sk 

f M 2 , and the definition r] k = 
= 2nnT with n = 0,1,2- ■■ is 
the meson Matsubara frequency in the imaginary time 



Zsk H 



k 2 ) 



■Ml 
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formalism of finite temperature field theory, and we have 
considered different rcnormalization constants Zg k and 
Zpk for the scalar and pseudoscalar fields. In order to 
complete the set of flow equations, we need equations for 
the evolution of Zsk and Zpfc|2(j. To the one- loop level 
they read 

d ^ z sk = ^§^E Tr (D 2 Sk r sss D 2 Sk r sss ) (20) 



6tt 2 

z pl k4 
6n 2 



spp , 



-l 

'Sk 



■M% 



Different from the 0(4) 
2 and V)et{(j)) + Det{(^) 



and the corresponding propagators V < 
and V p \ = u 2 n + k 2 + A 
model[||], while Tr((fS (f>), Tr( 
are invariant under the SUl(3) x SU p(3) transformation, 
not all the meson masses can be expressed with the eigen- 
values of these operator composites. 

We now expand the potential around the minimum 
a uk and a sk . Shifting the field a u = a uk + 5a uk and 
a s = a sk + 5a s k, the derivative of the potential on the 
left hand side of the flow equation could be expressed as 
powers of 5a uk and Sa sk , 



-d k Z p l 



z sk Z Pkk 4 
3tt 2 



where (r sss )^ c = d 3 U k /da a da b da c , (T S pp) bc = 
d 3 U k /da a dir b dir c and (Tfsp)^ = d 3 U k /dTr a da b dir c are 
9x9 matrices for the three-line vertexes with a fixed 
external meson a. In the above discussion we have as- 
sumed that the wave function renormalization constant 
depends only on the Lorentz transformation property of 
the mesons, but is independent of the detailed meson 
types. In the following numerical calculations related 
to the three-line vertexes, we take the Goldstone mode 
a = 4. 

As we will see from the numerical calculation in the 
next section, the contribution from the wave function 
rcnormalization to the thermodynamics of the system is 
very small and can be safely neglected as a first order 
approximation. In this case, by taking Zg k = Z Pk = 1, 
the flow equation (Til?]) is further simplified as 



d k U k 



6tt 2 



Tj2^(Dsk + D Pk ) 



(21) 



Note that the parameters h u and h s (or ho and hg) 
which explicitly break the chiral symmetry are scale in- 
dependent, and only the condensates a u and cr s (or ao 
and (jg) which spontaneously break the chiral symmetry 
depend on the scale k. In the treatment of the flow equa- 
tion with classical potential, the mass and coupling con- 
stant renormalization leads to four fc-dependent param- 
eters mic, Cfe, Aifc and \2 k controlled by the flow equation 
(|2ip , and the condensates a uk and a sk are determined by 
the gap equations (fTTI) . 

Considering the partially degenerated diagonal ele- 
ments and the simple off-diagonal structure of the mass 
matrices Ms and M P , the trace Tr in (f2~Tj) can be easily 
done, and the flow equation becomes 



d k U k = ^T^[Tr(X> Sfe +Cp fe ) 

n 

4D Kk ] 



67T 2 

n 

+ 3D aok +4D Kk + 3D nk 
with two 2x2 mixed matrices 



(22) 



Ms 
M P 



(M 2 ) 00 (AfJ)o8 



(Ms) 



s )so 



(M 2 



S )8S 



(M 2 )oo 
(M 2 P ) m 



(M 2 ) 08 

(M 2 )s8 



(23) 



= U k 



U k (a uk + 5a uk ,a sk + Sa sk ) 
dU k 



dSa m 



-5cr uk 



1 d 2 ii k 2 



f 



dU k 1 d 2 U k 2 

tt? oa sk + - 2 da uk 

d 2 u k 



2 dSa uk d5a sk 



Sa uk Sa 



sk 



(24) 



with the definition U k = d k U k . By comparing the co- 
efficients of Sa ukl 5<J sk , Sa^ k , Sa 2 k and 5a uk Sa sk on the 
left and right hand sides of the flow equation (f2"2"j) . one 
obtains four independent differential equations for the 
six parameters m 2 , c k , Xik, \2 k , o uk and a sk . Together 
with the two gap equations, their fe-dependence are fully 
determined. Note that in deriving the four flow equa- 
tions we have used the relations d k 8a uk = —d k a uk and 
d k 5a sk = —d k a sk . 

Before numerically solving the flow equations which 
will be done in the next section, we firstly discuss their 
chiral limit analytically. From the relation between /i 8 
and tig, 



ig = m K a s , 



(25) 



since there is no reason for the k meson to be massless in 
chiral limit, the condensate a% should vanish in the case 
of ho = /is = 0. In this limit, there remains only one gap 
equation for the condensate ao, 



c _ /. A 2 \ 

7r o+ l Ai+ T» <r ° 



0, 



(26) 



and the off-diagonal elements in the mass matrices Ms.p 
disappear. In the phase of spontaneous chiral symmetry 
breaking, there are eight pseudoscalar Goldstone modes 
7r, K and rj which dominate the thermodynamics of the 
system. From the flow equation (|22[) . the heavy modes 
do not contribute much to the flow, we can keep only 
the terms with the Goldstone modes. In this approxi- 
mation the flow equation for the chiral condensate oofc is 
decoupled from the others, 



dkcr 2 k 



8k 4 

2k 
3^ 



i 



2n k 4 



hfc 2 ) 2 

2^n k + 2^n 2 /T) (27) 



with the Fcrmi-Dirac distribution n k = 1/ (e 



k/T 



!)• 
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It is easy to see that the right hand side of the flow 
equation (|27p is always positive, which leads to a mono- 
tonically increasing aok with the momentum scale k at 
any temperature. Therefore, when we start at an ultra- 
violet momentum A, the physical condensate at k = 
is guaranteed to be finite from the evolution of the flow 
equation. On the other hand, from the meson frequencies 
u> n = 2mtT, the condensate drops down with increasing 
temperature at any momentum scale k, which may lead 
to a phase transition of chiral symmetry restoration at a 
critical temperature T c . In fact, the flow equation can be 
analytically integrated out with the solution 



'Oft 



(T) 



'OA 



(T) + f(A,T)-f(k,T) 



(28) 



with the definition 



k 



5|£ 

T 



12-ln (-n fc ) - 12Li 2 | . 



hi. 



k/T 



(29) 



remarkably by the temperature, the temperature depen- 
dence of the parameters at the ultraviolet momentum 
can be safely neglected. Therefore, we take the tempera- 
ture independent initial values m A (T) = m A (0),c A (T) = 
c a (0),Aia(T) = Ai A (0),A 2A (T) = A 2A (0),a„ A (T) = 
^ma(O) and (7 s a(T) = ct sA (0), and they are so chosen to 
reproduce their vacuum values at k = discussed in Sec- 
tion [2J by solving the flow equations and gap equations 
at zero temperature. 

FigJT] shows the evolution of the condensates a u k , <J s k 
and ask at zero temperature. With increasing momen- 
tum scale fc, the condensate a^k drops down continuously 
and approaches to zero at k ~ I GcV. From the relation 
hg = ra 2 Kk ask, the fc-independence of the chiral breaking 
parameter h% leads to a divergent meson mass m K k at 
k ~ 1 GeV. Therefore, we take 1 GeV as the ultravio- 
let momentum A for the evolution of the renormalization 
parameters. We have checked that the numerical results 
are insensitive to the value of the momentum cut A. Both 
the light and strange quark condensates a u k and a s k go 
up with increasing scale k, as we analyzed in the end of 
Section [5] in chiral limit. 



where Li 2 (a;) = Y^tLi xl /l 2 is the polylogarithm function. 
At very high momentum, the temperature effect on the 
system becomes weak, and we can reasonably take the 
boundary condition of the flow equation (Toa (T) at finite 
temperature as the one ctoa(O) at zero temperature. The 
condensate ctoa (0) is determined by reproducing the vac- 
uum value <7oo(0), 



'OA 



(0) 



'oil 



(0) + /(0,0)-/(A,0). 



(30) 



The critical temperature T c for the chiral phase transition 
is then determined by 



OA 



(0) + /(A,T c )-/(0,T c ) = 0. 



(31) 



In the symmetry restoration phase with T > T c , the 
potential expansion should be around the zero conden- 
sate (To = 0. In this case, all the mesons become degen- 
erate with mass m 2 and their contributions to the flow 
equation are the same. The procedure for the SU(2) 
model is discussed in 129, \2m ■ 



3. NUMERICAL RESULTS 

In this section we show our numerical results for 
chiral symmetry and J/a(1) symmetry restoration at 
finite temperature. In the case with explicit chi- 
ral symmetry breaking, one has to solve the coupled 
four flow equations, if the wave function renormaliza- 
tion is neglected, together with the two gap equa- 
tions. The initial condition for the four differential 
equations at a fixed temperature is the values of the 
six parameters at the ultraviolet momentum A, namely 
77i a (T),ca(T),A 1 a(T),A 2A (T), ( 7 u a(T) anda sA (T). Con- 
sidering the fact that the system at high enough mo- 
mentum is dominated by the dynamics and not affected 
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FIG. 1: (Color online) The evolution of the condensates 
&uki<Tsk an d o"8fc at T = 0, scaled by their vacuum values 
at k = 0. 

The temperature dependence of the light and strange 
quark condensates at k = is shown in FigJ^l While the 
both condensates are almost constants at low tempera- 
ture T < 50 MeV, they monotonically decrease at high 
enough temperature, as we expected from the analysis 
in chiral limit. Since strange quarks are much heavier 
than light quarks, which is reflected in the values of the 
explicit symmetry breaking parameters h u and h s , the 
SU(2) symmetry restoration should be much faster than 
the SU(3) symmetry restoration. This is the reason why 
the light quark condensate drops more rapidly than the 
strange quark condensate. The qualitative behave of the 
condensates shown here is similar to what obtained in 



the frame of CJT approximation 17 1. 

What is the effect of the wave function renormaliza- 
tion on the condensates? Solving the flow equations 
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50 100 150 200 250 

T (MeV) 

FIG. 2: (Color online) The temperature dependence of the 
light and strange quark condensates a u o, &s0, scaled by their 
vacuum values at T = 0. 



(|24[) for the renormalized mass and coupling constants 
mf , Cfe, X\k and \2k and ([20)) for the renormalized wave 
functions Zsk and Zpk and the gap equations f) 1 1 [) for the 
condensates a u k and <7 s fc, and taking the vacuum values 
Zso = Zpo = 1 at T = k = 0, the momentum scale de- 
pendence of Zg k and at zero temperature are shown 
in FigJ3J While the scalar part varies strongly with the 
scale k, the renormalization for the pseudoscalar mesons, 
among which the Goldstone modes dominate the ther- 
modynamics of the system, is a very smooth function in 
the whole region < k < A. Considering the fact of 
rjk oc duZ^ 1 in the flow equation (|19|) . the contribution 
from the wave function renormalization to the averaged 
action is expected to be slight. From FigJH the temper- 
ature dependence of the light and strange quark conden- 
sates, including the contribution from the wave function 
renormalization, is almost the same as the one shown in 
FigEI Therefore, for the following numerical calculations 
we will omit the wave function renormalization as a first 
order approximation. 

The Ua(X) anomaly is described by the topological 
susceptibility [|| 

X = J d 4 x(0\T(Q(x)Q(0))\0) (32) 

with Q(x) being the topological charge determined by 

the QCD coupling constant g and the gauge field tensor 
pa 

x fJ,V> 



Q(x) = ^e^F^(x)F^(x), (33) 

where t^ wpa is the antisymmetric tensor. By introducing 
the anomaly, the ninth Goldstone mode rj 1 is no longer 
massless. The topological susceptibility can be related to 
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0.8- 



0.4- 




200 400 600 800 1000 

k (MeV) 

FIG. 3: (Color online) The momentum scale dependence of 
the wave function renormalization constants at zero temper- 
ature. 
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FIG. 4: (Color online) The temperature dependence of the 
light and strange quark condensates, including the wave func- 
tion renormalization. 

the pseudoscalar meson masses in the linear sigma model, 




It is determined not only by the three-line coupling c 
but also by the four- line coupling A2, although the Ua (1) 
symmetry is broken only by the determinant term which 
is irrelevant to A2. 

The temperature dependence of the topological suscep- 
tibility is shown in FigJS] It decreases with temperature 
in the region of T > 100 MeV, indicating a continuous 
restoration of ?7a(1) symmetry. Different from the lattice 
simulation for the SU(3) Yang-Mills theory where the 
susceptibility drops down rapidly, the susceptibility re- 
mains still 30% of its vacuum value at high temperature 
T = 250 MeV. Unlike the condensates which never in- 
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crease with temperature, the susceptibility shows a slight 
increase in the low temperature region of T < 100 MeV, 
which is observed also in simulations in the interacting 
instanton liquid model0. 



1.2 




0.2 , 1 , 1 . 1 , 1 , r- 

50 100 150 200 250 
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FIG. 5: (Color online) The temperature dependence of the 
topological susceptibility %■> scaled by its vacuum value at 
T = 0. 

The chiral symmetry and ?7a(1) symmetry restoration 
shown above control the meson masses in hot medium. 
The temperature dependence of the meson masses is 
shown in Figs [6] and [7l There are four different mesons for 
each species, the triplet ao, quartet k, and the mixed a 
and /o for the scalar mesons, and the triplet 7r, quartet K, 
and the mixed rj and rj' for the pseudoscalar mesons. As 
temperature increases, the scalar meson masses m/ , m K 
and m ao drop down monotonically, but /o and k are heav- 
ier than ao at T < 200 MeV and become lighter than <z 
at T > 200 MeV. The Goldstone modes ir and their chiral 
partner a become degenerate at high temperature, due 
to the chiral symmetry restoration. The similar behave 
is for the mesons ao and rf . At T = 250 MeV there 
are ~ m a ~ 200 MeV and m ao ~ m v > ~ 900 MeV. 
The mixed mode // behaves similarly to the other seven 
Goldstone modes ir and K. While the mass of the other 
mixed meson rj slightly goes up at low temperature and 
drops down at high temperature, the difference between 
r\ and rj decreases with increasing temperature, due to 
the partial restoration of the ?7a(1) symmetry 

In chiral limit, the condensate a& disappears and the 
flow equation for the condensate (To is decoupled from 
the others and can be analytically solved as (|28|). The 
evolutions of CTofc at different temperature is shown in 
FigJHl While the temperature dependence at low mo- 
mentum scale k is remarkable, the condensate is almost 
T-indcpendent when the scale is large enough. This sup- 
ports our choice of T-indcpcndcnt initial condition at 
k = A. With increasing temperature, the solution at 
k = approaches to zero continuously, indicating the 
phase transition of chiral symmetry restoration. The 
temperature dependence of the physical condensate and 
topological susceptibility at k — is shown in FigjS] 




50 100 150 200 250 

T (MeV) 



FIG. 6: (Color online) The temperature dependence of the 
scalar meson masses. 
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FIG. 7: (Color online) The temperature dependence of the 
pseudoscalar meson masses. 

The condensate drops down much faster than that in 
the general case with explicit chiral symmetry breaking, 
see Figl^J The critical temperature for the chiral phase 
transition determined by aoo(T c ) — is T c = 130 MeV. 

The topological susceptibility (IMl) for the JZa(1) 
anomaly is simplified as 

X - ^cal (35) 

in chiral limit. After a rapid decrease in the vicinity of 
vacuum, it becomes very smooth and finally vanishes at 
the same critical temperature for the SUa(3) symmetry 
restoration. 



4. CONCLUSION 

We investigated in this paper the chiral symmetry and 
[/^(l) symmetry restoration at finite temperature, by ap- 
plying the functional rcnormalization group to the SU (3) 
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FIG. 8: The evolutions of the condensate oofc at different 
temperature in chiral limit. 
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linear sigma model. Wc derived the flow equations for 
the mass, coupling and wave function renormalization 
parameters in the local potential approximation and the 
two gap equations for the light and strange quark con- 
densates. 



In chiral limit, wc analytically solved the decoupled 
flow equation for the chiral condensate and analyzed 
its momentum scale and temperature dependence. The 
eight Goldstone modes are guaranteed in vacuum and at 
finite temperature before the chiral restoration, and the 
two phase transitions for the chiral symmetry and Ua{1) 
symmetry restoration take place at the same critical tem- 
perature T c = 130 MeV. 



In general case with explicit chiral symmetry breaking, 
we numerically solved the coupled flow and gap equa- 
tions at finite temperature, starting from the classical 
potential at the ultraviolet momentum k = A ~ 1 GeV 
and extracting physics including quantum fluctuations at 
k = 0. In this case a partial restoration of the SUa(3) 
and C/a(1) symmetries is observed. Different from the re- 
sults obtained in other approximations like CJT method, 
the light and strange quark condensates drop down with 
temperature slowly. As a result of the partial restoration 
of the two symmetries, the pseudoscalar triplet 7r and its 
chiral partner a ( the scalar triplet ao and 77') become 
degenerate at high temperature, and the difference be- 
tween the mixed modes r\ and rf gradually disappears in 
the limit of high temperature. 



FIG. 9: The temperature dependence of the chiral condensate 
<To and topological susceptibility x m chiral limit, scaled by 
their vacuum values. 
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